o 
o 
o 

(N 
> 

o 






Monoclinic and triclinic phases in higher-order Devonshire theory 

David Vandcrbilt and Morrel H. Cohen 
Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08854-8019 

(November 16, 2000) 

Devonshire theory provides a successful phenomenological description of many cubic perovskite 
ferroelectrics such as BaTiOs via a sixth-order expansion of the free energy in the polar order 
parameter. However, the recent discovery of a novel monoclinic ferroelectric phase in the PZT 
system by Noheda et al. (Appl. Phys. Lett. 74, 2059 (1999)) poses a challenge to this theory. Here, 
we confirm that the sixth-order Devonshire theory cannot support a monoclinic phase, and consider 
extensions of the theory to higher orders. We show that an eighth-order theory allows for three kinds 
of equilibrium phases in which the polarization is confined not to a symmetry axis but to a symmetry 
plane. One of these phases provides a natural description of the newly observed monoclinic phase. 
Moreover, the theory makes testable predictions about the nature of the phase boundaries between 
monoclinic, tetragonal, and rhombohedral phases. A ferroelectric phase of the lowest (triclinic) 
symmetry type, in which the polarization is not constrained by symmetry, does not emerge until 
the Devonshire theory is carried to twelfth order. A topological analysis of the critical points of the 
free-energy surface facilitates the discussion of the phase transition sequences. 

PACS: 77.80.Bh, 77.84.Dy, 64.70.Kb 
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I. INTRODUCTION 

Many of the most important and best studied ferroelec- 
tric materials adopt the simple-cubic perovskite struc- 
ture at high temperature and undergo structural phase 
transitions to distorted ferroelectric structures at lower 
temperature. Among the best known simple compounds 
of this kind are BaTiOs and PbTiC>3. Upon cooling, 
BaTiOa undergoes a sequence of ferroelectric transitions, 
first from the cubic (C) to a tetragonal (T), then to an 
orthorhombic (O), and finally to a rhombohedral (R) 
phase. Passing through this sequence, the polarization 
P first vanishes in the C phase, and then becomes ori- 
ented in the [001], [Oil], and [111] directions in the T, O, 
and R phases, respectively. PbTiC>3 undergoes a single 
ferroelectric transition from the C to T. phase. 

In a classic 1948 paper, Devonshire^ was able to ex- 
plain the observed phases and phase transition sequence 
quite naturally in terms of a phenomenological Landau- 
type expansion of the free energy in terms of the ferro- 
electric order parameter P. Making use of cubic sym- 
metry and truncating the expansion to sixth order in P, 
Devonshire was able to arrive at a simple model with 
only a single temperature-dependent second-order coeffi- 
cient, and only three temperature-independent higher- 
order coefficients. The polarization P is the primary 
order parameter, and the crystallographic labels (T, R, 
etc.) refer to the distortions induced by the polariza- 
tion and the resulting strain. Despite its simplicity, this 
model could successfully reproduce the phase-transition 
sequenceEl and the piezoelectric and other properticsa of 
BaTiC>3. With a simple modification of the anharmonic 
coefficients, the qualitative behavior of PbTiC>3 could be 
equally well reproduced. 

However, the material that is currently in most wide- 
spread use for piezoelectric transducer and related ap- 



plications is the solid solution PbZri-^T^Os, commonly 
known as PZT. The standard understanding of the phase 
diagram of PZT has been as followsfl PZT undergoes a 
transition from the simple cubic C phase to a ferroelectric 
phase at a Curie temperature T c that ranges from about 
490°C at x=l to 230°C at x=0. The transition occurs 
to the T phase for x greater than about 0.48, and to the 
R phase for smaller x. At T < 150°C or x < 0.1, some 
more complex phases involving antiferroelectric (AFE) 
or antiferrodistortive (AFD, i.e., involving rotations of 
oxygen octahedra) displacements also occur. The phase 
boundary between the simple T and R phases, known as 
the morphotropic phase boundary (MPB), is an almost 
vertical line in the x-T plane at about x=0.48. Haun et 
al.a have successfully extended the Devonshire model to 
the case of PZT by including the AFE and AFD degrees 
of freedom in the phenomenological free energy, still in- 
cluding only terms up to sixth order overall. This model 
successfully described the simple R-T transition that was 
understood to occur below T c at the MPB. 

As it happened, a surprise was in store. By working 
with highly purified and carefully prepared samples of 
PZT, Noheda and coworkers have recently shownETEl that 
a sliver of monoclinic (M) phase actually interposes itself 
between the R and T phases in a very narrow composi- 
tion range (of order 3-4% in x). That is, at least below 
^100°C, the transition is first from T to M at an x C 2 be- 
tween about 0.48 and 0.51 (depending on temperature), 
and then from M to R at x c \ ~ 0.47, with decreasing x. 
The orientation of P is, respectively, along [001], [uuv] 
(u < v), and [111] in the T, M, and R phases, respec- 
tively. The experiments have not yet clarified whether 
or not a direct T-R transition occurs in the higher tem- 
perature range 100°C< T < T c , or whether the sliver 
of M phase instead persists up to the Curie tempera- 
ture T c ~ 370°C. Neither the Devonshire theoryEl nor 



the modification of Haun et alp predicted the possible 
occurrence of the M phase. On the other hand, simu- 
lations based on a first-principles effective Hamiltonian 
approachE] have very recently provided confirmation of 
the existence of_the M phase in just the observed com- 
position range.E£l 

We make a brief aside to establish notation. In the 
C phase (space group Pm3m), P = 0. When P is con- 
strained to a symmetry axis lying along [001], [111], or 
[Oil], the resulting phase and space-group labels become 
T (P4mm), R (i?,3m), or O (Amm2), respectively. Sim- 
ilarly, the M phases arise when P is confined to a mir- 
ror plane. We can distinguish three cases: Mq (Pm), 
in which P is along [Ouv]; and Ma and Mb (both Bm, 
sometimes also denoted Cm) in which P lies along [hud], 
with u < v and u > v, respectively. The newly observed 
phase of PZT is of type Ma-lHj If P is unconstrained by 
symmetry, strain coupling leads to a triclinic phase (PI). 
This exhausts the possible reduced-symmetry states of a 
cubic perovskite crystal generated by the emergence a 
single non-zero ferroelectric order parameter, although 
other phases can occur if AFE and AFD distortions are 
also present. 

Relatively few examples are known of low-symmetry 
ferroelectrics in which P is only constrained to a symme- 
try plane, or in which P is unconstrained by symmetry. 
The discovery of the Ma phase in a cubic perovskite is 
thus of considerable note, even aside from the fact that 
it had been missed for so long and aside from the poten- 
tial importance of this phase far understanding the large 
piezoelectric response in PZT.fl 

The failure of the phenomenological theories of Devon- 
shircfl and Haun et aln to describe the existence of the 
observed Ma phase raises an interesting question: What 
is the simplest and most natural phenomcnological model 
that does predict such, a monoclinic phase? Recent work 
of Souza Filho et alEA confirms that the sixth-order De- 
vonshire expansion does not allow for the occurrence of a 
monoclinic phase and suggests instead a model in which 
the vanishing of a shear elastic constant at the critical 
temperature drives the transition to the Ma phase. How- 
ever, this requires the introduction of an additional in- 
stability that is unrelated to the ferroelectric oneJi3 That 
this is unnecessary is demonstrated by the work of Bel- 
laiche, Garcia and Vanderbilt.EHl In their simulations, the 
shear moduli are taken as temperature- independent, and 
yet the transition to the monoclinic phase still occurs, 
driven by the tilting of the ferroelectric polarization away 
from the symmetry axis. 

Suspecting instead that the problem is simply related 
to the truncation of the expansions of Devonshirea and 
Haun et ala to sixth order, we consider the addition of 
terms of eighth and higher order to the Devonshire model. 
The following questions then arise. At what order in the 
expansion do monoclinic phases first appear in the phase 
diagram? And, for that matter, at what order do triclinic 
phases first appear? 

The phase diagram of such a model consists of fields 



in the parameter space of the model (labeled C, T, O, 
R, M, etc.) within which the order parameter has the 
specified symmetry (that is, the absolute minimum of 
the free energy occurs at a point in the order-parameter 
space of that symmetry). On a trajectory in the phase 
space which crosses a phase boundary, that global min- 
imum can change into a local minimum, a saddle point, 
or a maximum, or simply disappear. Understanding how 
the set of critical or stationary points of the free energy 
(its minima, saddle points, and maxima in the order- 
parameter space) varies with the parameters of the model 
can thus be an important aid in understanding the phase 
diagram. So, in addition to searches for the global min- 
ima by direct computation, we use here a topological 
analysis of the Devonshire theory to elucidate the critical- 
point structure of the free energy and answer the above 
questions. 

We show here that the simplest extension of the Devon- 
shire theory, to eighth order in the ferroelectric order pa- 
rameter, naturally admits all three kinds (Ma, Mb, Mq) 
of monoclinic phases. Moreover, the eighth-order model 
makes specific and testable predictions for the types of 
phase boundaries that can occur. For example, it pre- 
dicts that the R-Ma and T-Ma transitions should be of 
first and second order, respectively. On the other hand, 
we find that the model has to be extended all the way to 
twelfth order in order to describe a triclinic ferroelectric 
phase. 

We therefore suggest that the most natural explanation 
for the occurrence of the Ma phase in PZT is that the 
free-energy surface is unusually anharmonic in this mate- 
rial, such that the eighth-order terms play an important 
role. It is not difficult to speculate why this might be 
the case. First, PZT is a disordered material; averag- 
ing over the chemical disorder, which plays the role of 
a quenched random field, may tend to generate higher 
orders in the phenomenological energy expansion. Sec- 
ond, as will be discussed in Sec. M, there is considerable 
evidence that the ferroelectric transitions in PZT have 
a strong order-disorder character; mapping onto a dis- 
placive picture may then also tend to generate higher- 
order terms. And finally, of course, thermal fluctuations 
and coupling to strain may play some role. In any case, 
we show below that once one accepts this simple hypoth- 
esis of the importance of eighth-order terms in the free- 
energy expansion, then the hehavior of the monaeJinic 
phase observed in experimentsBlffl and simulations!^ can 
be understood quite naturally. 

This paper is organized as follows. Section O estab- 
lishes the notation used for the expansion of the free en- 
ergy, and reviews the symmetry considerations that lead 
to a simplified form of this expansion. The rules gov- 
erning the numbers and types of stationary points that 
may occur in the order-parameter space are reviewed in 
Sec. III. The behavior of the models obtained by trun- 



cating the expansion at higher and higher order (fourth, 
sixth, eighth, tenth, and twelfth) are then carefully elu- 
cidated in Sec. IV. Section [V| gives a brief discussion of 



a microscopic model in which the T and R phases can be 
regarded as arising from fluctuations among neighbor- 
ing, symmetry-equivalent local Ma states. Finally, we 
conclude with a brief summary and discussion of future 
prospects in Sec. [VJ. 

II. FORMALISM 

We consider the case of a structural phase transition 
governed by a single continuous vector order parameter 
u, such that the free energy 



F(u,<T i ,T) = E(u,r H ,S)-J2°iVi-TS 



(1) 



is symmetric with respect to operations of the cubic point 
group. Here o~i and r)i are the stress and strain tensors 
in Voigt notation, and T and S are temperature and en- 
tropy, respectively. We have in mind primarily the case 
in which u is the ferroelectric polarization P in a mem- 
ber of the cubic perovskite class, such as BaTiC>3 or PZT, 
but the use of the symbol u serves as a reminder that the 
formalism applies to a variety of other cases as well. For 
a crystal with stress-free boundary conditions at some 
given temperature, we can expand 

F(u)= Cooo + C* 20 o (ul +u 2 y + u 2 z ) + C 400 (u* + u A y + uj) 
+C220 {u 2 x u 2 y + u\u\ + u 2 y ul) + C 60 o (ul + m^ + ul) 



+C420 {u x [u 2 y H 
+C222 u 2 x u 2 y u 2 z 



+ u y [u 2 x +u 2 z ] 



■ ut \ui 



(2) 



where terms up to sixth order in u have been written 
explicitly, and the coefficients have been renormalized to 
subsume the couplings to strain. 

For the ferroelectric phases we simplify matters further 
by focusing on the energy as a function of the orienta- 
tion of the vector order parameter and introducing the 
function 



G(u) 



min F(\i) 

u||u 



(3) 



Thus, G(u) represents the ground-state energy subject to 
the constraint that the order parameter has given orien- 
tation. In most cases of interest, it is reasonable to expect 
that G(u) will be a smooth function of u£j In this case, 
and suppressing the uninteresting constant term in the 
expansion, it follows that 



G(u) 



a 4 54-|-a 6 5 6 4 
+ aio5456 



+- 01204 



Oi 2 5 6 + 



where 



54 (u) = x 2 y 2 + x 2 z 2 + y 2 z 2 
5 6 (u) = x 2 y 2 z 2 



(4) 



(5) 



TABLE I. Values of the cubic invariant functions defined in 
Eqs. (Eh and (pi), evaluated at tetragonal (T), rhombohedral 
(R), and orthorhombic (O) orientations of the order parame- 
ter u. 
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and u = (x, y, z) with x 2 
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1. All independent, 



symmetry-allowed terms up to twelfth order are explicitly 
given in Eq. (0); higher orders will not be needed here. 

The particular form of the cubic invariants appearing 
in Eq. (3) is largely arbitrary. For example , one could use 
the "kubic harmonics"li-3 instead. In Sec. IV D , we shall 
make use of the expansion of Eqs. (0-0) above. However, 
for the presentation of numerical results to be given in 



Sec. 1VC 
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A(u) 


= 1254 






U(u) 


= 454- 


365 6 




/s(u) 


= 48 54 2 


-1254 



cs/s 



(6) 



365 6 



(7) 



Aside from normalization constants, which are chosen to 
make the range of each function roughly of order unity, 
this choice can be uniquely defined by the following re- 
quirements, (i) The invariant /„ contains no terms of or- 
der higher than n. (ii) All three invariants vanish identi- 
cally for a "tetragonal" (T) value of the order parameter, 
e.g., u = (100). (hi) f& and /§ still vanish for a "rhombo- 
hedral" (R) value of the order parameter, e.g., (lll)/\/3. 
(iv) /§ alone vanishes for an "orthorhombic" (O) value of 
the order parameter, e.g., (HO)/-^. The values of these 
functions evaluated at symmetry directions are summa- 
rized in Table I for later reference. Requirement (ii) just 
reflects an arbitrary choice of zero; requirements (iii) and 
(iv) simplify some later discussion. For example, phase 
boundaries at which G(T), G(R), and G(O) become de- 
generate in pairs are easily located (see Sec. IVC). 

A few comments about the transition from -F(u) to 
G(u) in Eq. ([|) are in order. First, our "reduced" the- 
ory in terms of G(u) can only describe the transitions 
among the ferroelectric phases, not the transition to the 
high-temperature cubic phase (u = 0). Second, the min- 
imization over the magnitude of u in Eq. (H) generally 
introduces higher orders into the reduced theory. For ex- 
ample, truncation of Eq. (0) at sixth order may still lead 
to terms of eighth and higher order in Eqs. (Q) and (g). 
Thus, a statement about what kinds of critical points 
can occur simultaneously in an n th -order version of the 



reduced theory does not necessarily carry over to the n th - 
order version of the standard Devonshire theory. Nev- 
ertheless, a statement that a certain minimum order of 
expansion is needed for the existence of a certain kind of 
stationary point in the reduced theory does carry over to 
the standard Devonshire theoryEj For example, we show 
later that the reduced theory must be carried to eighth 
order to allow for a monoclinic minimum, i.e., an equilib- 
rium u lying in a (110) plane. Then the same statement 
applies to the standard Devonshire theory. For, suppose 
that Uq is the magnitude of u at the minimum, and let 
£?o(u) = F(uqu); then Go must have a monoclinic mini- 
mum. But Go is also of no higher than eighth order.-and 
so admits only the same type of minima as does G.I13 



III. TOPOLOGICAL CONSIDERATIONS 

The rules governing the numbers and types of criti- 
cal or stationary points of a scalar function defined on a 
continuous manifold emerge from a branch of algebraic 
topology called Morse theory|-ar, equivalently, the calcu- 
lus of variations in the largeJla These rules-have proved 
useful in analyzing lattice vibration spectraliU and the ex- 
citation spectra of crystalline materials generally. They 
have not, however, enjoyed the broad application in con- 
densed matter physics that might have been expected 
from these early successes. 

The rules take the form of a set of inequalities and one 
equality. The latter is the most powerful, and we confine 
ourselves to considering only it explicitly. G(u) is an 
analytic function defined on the surface of a sphere, a 
closed two-dimensional manifold of genus zero, for which 
the equality becomes 



(a) 



N - JVi + N 2 = 2 



(8) 



where Nq, N±, and N2 are the number of minima, saddle 
points, and maxima, respectively. This formula assumes 
only analyticity in G(u) and that its Hessian has no van- 
ishing eigenvalues at the critical points. 

Because the function G(u) has cubic symmetry, any 
stationary points that occur must be members of 
symmetry-related families of stationary points. Thus, it 
is natural to focus attention on an irreducible wedge cor- 
responding to 1/48 of the unit sphere, in terms of which 
Eq. (pt) can be rewritten as 



E' 



ili 



(9) 



where j runs over the stationary points located in the in- 
terior or on the boundary of the irreducible wedge, Uj is 
a degeneracy factor counting the number of images of the 
stationary point generated by the cubic symmetry group 
operations, and jj is +1 for a maximum or minimum and 
— 1 for a saddle point. Symmetry requires only that there 
be stationary points at T, O, and R; those at T and R 





FIG. 1. (a) Contours of constant G(u) = — /g(u) on the 
unit sphere, (b) Labels for symmetry points and lines of the 
irreducible wedge, (c) Corresponding degeneracy factors n 
indicating number of images on the full unit sphere, (d) Sum- 
mary representation of the behavior of G(u), in which open, 
shaded, and filled symbols represent maxima, saddle points, 
and minima of G, respectively. Arrows indicate "downhill" 
flow lines. 



must be maxima or minima while that at O can be of any 
type. Such a symmetry set of stationary points may or 
may not be large enough to satisfy the Morse relations as 
well. A set which contains the smallest number of criti- 
cal points satisfying both the topological and symmetry 
requirements is denoted a minimal set .113 

This type of analysis is illustrated in Fig. 1. Figure 1(a) 
shows an arbitrarily chosen function G(u) corresponding 
to C4 = eg = and cs = — 1 in Eq. (ra), having maxima 
at T, R, and O; saddle points on the lines connecting T- 
O and O-R; and a minimum on the segment connecting 
T-R. Figure 1(b) specifies the notation that we shall use 
to identify the symmetry points and lines. These are 
chosen to correspond to the labels of distorted crystal 
structures (tetragonal, rhombohedral, or orthorhombic 
for the symmetry points; monoclinic of type 'A,' 'B,' or 
'C for the symmetry lines; and triclinic for the case of u 
pointing to the interior of the wedge). The corresponding 
degeneracy factors n are given in Fig. 1(c). The overall 
behavior of the function inside the wedge is summarized 
in Fig. 1(d), in which the stationary points are marked 
by symbols that are open, shaded, or filled for maxima, 



(a) T 




(a) M 




(d) R 




(d) Tri 



FIG. 2. Sample topologies giving rise to tetragonal (a-c), 
rhombohedral (c-e), and orthorhombic (e-f) phases. Open, 
shaded, and filled symbols represent maxima, saddle points, 
and minima, respectively. 



saddle points, or minima, respectively. In this particular 
example, it can easily be seen that Eq. (|[) is satisfied 
(6+8+12-24-24+24=2). 

Using Eq. (|9|), one can begin enumerating the possi- 
ble topological diagrams for the cubic system. Figs. 2(a) 
and 2(d) show the two possible diagrams in which there 
are only three stationary points in the irreducible wedge, 
the only possible minimal sets. In Fig. 2(a) the function 
has a unique minimum at T, so the crystal ground state 
would be tetragonal. A similar situation holds for the 
rhombohedral case of Fig. 2(d). Figures 2(b-c) and (e-f) 
show four of the six configurations that can exist when 
there are exactly four stationary points in the wedge. 
(Two others, similar to 2(b) and 2(e) but with the the 
role reversals T<->R and Mb^Mc, are not shown.) In 
Fig. 2(c) there are two local minima, so that the ground 
state could be of type T or R depending on which min- 
imum is deeper; a first-order transition between T and 
R may occur by a crossing of the minima. The simplest 
configuration leading to an orthorhombic ground state is 
shown in Fig. 2(f). 

In order to arrive at configurations corresponding to 
monoclinic or triclinic phases, more than four stationary 
points are required. Figs. 3(a-c) illustrate possible con- 
figurations for monoclinic phases of type A, B, and C, 
in which the order parameter lies on the T-R, R-O, and 
O-T lines, respectively. Figure 3(d) illustrates a possible 
triclinic phase. It is straightforward to check that Eq. (Q) 
is satisfied for each configuration in Figs. 2 and 3. 

Of course, if the expansion of Eq. (|4j) or (|6j) is trun- 
cated at a certain order, the possible types of topological 
behavior will be limited by the enforced "smoothness" of 
the functions allowed at that order. The purpose of the 
following Section is to explore precisely this issue, i.e., 
to clarify what types of phases and phase transitions can 
occur at each order in the expansion. With this infor- 
mation in hand, one can then easily determine what is 
the minimal model needed to study a particular physical 
phenomenon of interest. 




FIG. 3. (a-d) Sample topologies giving rise to monoclinic 
phases A, B, and C, and the triclinic phase, respectively, (e) 
Degenerate minimum (heavy solid curve) that can occur when 
j3 — 37r/4. (f) Degenerate minimum (heavy solid line) that 
occurs for a = tan -1 (3) and f3 = n. Open, shaded, and 
filled symbols represent maxima, saddle points, and minima, 
respectively. 



IV. RESULTS 

A. Fourth-order theory 

If the expansion (H) is truncated at fourth order, the 
only non-constant cubic invariant is f± of Eq. (ffl). The 
T or R phase is favored as in Fig. 2(a) or 2(b) if C4 > 
or C4 < 0, respectively. A transition between T and R 
phases occurs at C4 = 0, but this transition is unphys- 
ical because the energy surface is perfectly flat at the 
transition. This degenerate behavior is an artifact of the 
truncation to fourth order. 



B. Sixth-order theory 

At sixth order in the expansion (|^), the behavior is 
governed by the two coefficients C4 and cq. Clearly a 
common rescaling of the coefficients by a positive scale 
factor is irrelevant, so without loss of generality we may 
set c| + Cg = 1 . It is then convenient to let 



C4 = cos(a) , cq — sin(a) 



(10) 



By doing numerical calculations and making plots such 
as that shown in Fig. 1(a), we find that only three phases 
exist at this level. For — tan -1 (3) < a < ir/2, the system 
is in the T ground state; for ir/2 < a < 57r/4, the system 
is in the R ground state; and for 57r/4 < a < 2ir — 
tan _1 (3) the system is in the O ground state. Fig. 2 
illustrates the sequence of topologies traversed as a is 
increased. Starting from a = — tan _1 (3), one finds T 
phases corresponding to Figs. 2(a-c) consecutively, until 
a first-order transition occurs to the R state at a = tt/2. 
Then one finds R phases as illustrated in Figs. 2(c-e) 
until a first-order transition occurs to the O phase at 



a = 57r/4. Finally, O phases corresponding to Figs. 2(e- 
f) are found as a is increased further up to 2ir — tan _1 (3). 
The transition from the O to the T phase is degenerate 
in the sense that the energy surface becomes exactly flat 
along the entire Mc symmetry line at the critical a (this 
being an artifact of truncation to low order). 

It is important to emphasize that that no monoclinic 
phase is possible in the sixth-order model. In the range 
7T + tan -1 (3/2) < a < 27T-tan- 1 (3) (i.e., 1.3128tt < 
a < 1.60247r) in which such a ID local minimum appears 
along the Ma symmetry line, it is always unstable in the 
second dimension (i.e., it is a saddle point), and the true 
minimum is at the O point as in Fig. 2(f). 

Further details of the sequence of topologies and the 
boundaries between them is given in the Appendix. 



C. Eighth-order theory 

When the model of Eq. (ph is carried to eighth order, it 
no longer becomes profitable to enumerate every possible 
topology, as was done above for the sixth-order theory. 
Instead, we choose to focus just on the "phase diagram" 
that is generated by finding the ground-state symmetry 
as a function of the parameters C4, eg, and eg. A common 
scaling of the magnitudes of these coefficients is again 
unimportant, so we can describe the phase diagram in 
terms of two dimensionless parameters that we may take 
as 

C4 = cos(a) , 

c 6 = sin(a) cos(/3) , 

c 8 = sin(a) sin(/3) , (11) 

where < a < it and < (3 < 2ir. From this point 
of view, the "parameter space" is just the unit sphere 
determined by polar and azimuthal angles a and j3 re- 
spectively. 

Figure 4 shows the phase diagram that emerges from 
a careful numerical study of the minimization of Eq. (0) 
as a function of a and (3. The plot is a mapping of the 
unit sphere onto the page. The points at the "north" and 
"south" poles (a — and a = it) are the only ones acces- 
sible in the fourth-order theory (Sec. IV A); the dotted 
vertical lines at (3 = and (3 = ir correspond to the locus 
of points in parameter space that were explored by the 
sixth-order model (Sec. IV B). 

As can be seen from Fig. 4, six of the seven possi- 
ble phases (i.e., possible symmetries of a non-zero order 
parameter) are accessed by the eighth-order model. In 
addition to the T, R, and O phases that appeared al- 
ready at sixth order, all three monoclinic phases (Ma, 
Mb, and Mc) are now stable in some region of the phase 
diagram. However, the areas covered by the Mc and es- 
pecially Mb regions are relatively small, so these phases 
may be harder to find in real systems than the Ma phase. 
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FIG. 4. Phase diagram in the space of parameters a and /3 
of the eighth-order theory as defined in Eq. ([ill). Solid and 
dashed lines are first-order and second-order phase bound- 
aries, respectively. The vertical grey line and the grey dots 
indicate cases for which degenerate minima occur. Fine verti- 
cal dotted lines indicate the domain of the sixth-order theory. 



Solid and dashed phase boundaries indicate transitions 
of first and second order, respectively, as determined nu- 
merically. According to Landau theory, there are two 
necessary conditions for a transition to be of second or- 
der: (i) the symmetry group of one phase must be a sub- 
set of the symmetry group of the other, Q C Qo; and (ii) 
from the order-parameter displacements that lead from 
the high-symmetry to the low-symmetry phase, it should, 
be impossible to construct a third-order invariant of £?o-" 
Transitions of type T-R, T-O, and O-R are necessarily 
first-order because of condition (i) , and transitions from 
the R phase to Ma or Mb are first-order because of (ii) .EJ 
As can be seen from the figure, these are precisely the 
boundaries that were found to be of first order. All oth- 
ers are found to be of second order, with the exceptions of 
the Ma-Mb and Ma-Mc boundaries, which form a ver- 
tical line at (3 — 37r/4 indicated by grey shading in the 
figure. Along this line, one finds a degenerate minimum 
connecting Ma and Mc phases as shown in Fig. 3(e) for 
0.2597T < a < 0.6407T, and a similar degenerate minimum 
connecting Ma and Mb phases for 0.6407T < a < 0.7407T. 
The degenerate behavior can be traced to the fact that 
ag = along the line (3 = 37r/4. The triple point connect- 
ing the T, O, and M c phases at a = tan" 1 (1/3) = 0.102tt 
and (3 = it is also a point at which a degenerate minimum 
occurs, as shown in Fig. 3(f). The degenerate behaviors 
are artifacts of the eighth-or der t runcation, as will be 
explained more fully in Sec. IVD. Finally, the reader 



is reminded that because the theory is based on a sin- 
gle polar orientational order parameter, Eq. (J3j), we do 
not have the ability to describe transitions to or from 
the paraelectric C phase, nor can we describe the more 



complex AFE or AFD phases. 

The variation of some physical variable, such as tem- 
perature, composition, or pressure, will correspond to a 
variation of the parameters a and (3 of the model in a 
way that is not easy to predict a priori. For BaTiC>3 
and KNDO3, for which the observed phase transition se- 
quence is R-O-T-C with increasing temperature, it must 
be the case that the system traverses a roughly vertical 
trajectory on the phase diagram of Fig. 4, somewhere 
in the range it < (3 < 3tt/2. (Insofar as eighth-order 
terms are small for these systems, the trajectory should 
be near (3 — it.) The R-O-T sequence can be visualized 
as traversing Figs. 2(d-f) and 2(a), in that order. 

On the other hand, systems such as PZT 
(PbZri_ x Ti x 03) that exhibit a morphotropic phase 
boundary (MPB), i.e., an R-T transition as a function 
of composition x, evidently cross the first-order R-T 
phase boundary in the vicinity of /3 = with increas- 
ing x. If this trajectory passes to the right of the triple 
point connecting R, T, and Ma phases at a = it/2, 
(3 = tan" 1 (1/3) = 0.102tt in Fig. 4, then the phase tran- 
sition sequencerJp«comes R-Ma-T, as recently observed 
experimentally.laCI'El The narrowness of the range of Ma 
phase, only a few percent in x, suggests that the tra- 
jectory passes rather close to the triple point. In fact, 
there are strong experimental indications of the possible 
existence of a triple point in the x-j^C jphase diagram of 
PZT near x ~ 0.47 and T ~ 100°Cmi Thus, it may be 
that the behavior near the triple point can be explored 
experimentally in the PZT system. 

Using the topological analysis introduced earlier, we 
can now clarify the nature of the T-Ma~R transition se- 
quence near the triple point. Referring to Fig. 5, we imag- 
ine traversing a downward trajectory of increasing a at 
fixed f3, slightly to the right of the triple point. Starting 
deep in the T phase, Fig. 2(a), a saddle point detaches 
from the O point and traverses the Mb line toward the 
R point, Fig. 2(b). After it passes through the R point 
and emerges on the "other side" (on the Ma line), we 
find ourselves in the situation of Fig. 5(a), the R point 
having been converted to a local minimum. Up to this 
point, the global minimum remains at T. Next, the T 
point converts from a local minimum to a local maxi- 
mum, with the simultaneous emission of a saddle point 
along Mc and a local minimum along Ma, as shown in 
Fig. 5(b). This event corresponds to the second-order 
T-Ma transition. The first-order Ma~R transition then 
occurs by the crossing of the energies of the local minima 
of Fig. 5(b). Once in the R phase, the local minimum 
and saddle point on the Ma line annihilate one another 
to give Fig. 5(c). Finally, deep in the R phase, the Mc 
saddle point eventually arrives at O, giving rise to the 
situation of Fig. 2(d). 

It should thus be emphasized that the present theory 
makes a definite prediction about the nature of the tran- 
sitions that occur in the T-Ma~R transition sequence. 
We can predict that, if it were possible to scan with de- 
creasing x at a temperature below that of the triple point, 
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FIG. 5. Topologies encountered in the T-Ma-R transition 
sequence near the triple point, (a) Topology of T phase near 
the second-order T-Ma transition, (b) Topology of Ma phase 
and, as well, of R phase near the first-order Ma-R transition, 
(c) Topology occurring deeper in the R phase. 



one would first find a continuous rotation of the polariza- 
tion from [001] into the (110) plane starting at a critical 
x c2 , and then a discontinuous jump to the [111] direction 
when there is a crossing of the free energies of the Ma and 
R phases at x c \ . Unfortunately, the fact that the T-Ma 
and especially the R-Ma boundaries lie almost vertically 
in the experimental x-T plane may make it difficult to 
test this prediction, since x can only be varied by prepa- 
ration of multiple samples. Nevertheless, this scenario 
seems to be supported by the numerical simulations of 
Ref. |l|. 

It is important to note that a triclinic phase does not 
occur anywhere in the phase diagram of the eighth-order 
model. However, we do note the possibility of observing 
new monoclinic phases of types Mb and Mc in a region 
near (3 ~ 0.87T. In fact, the sixth-order model (verti- 
cal dotted line at (3 = 7r) comes very close to yielding a 
monoclinic Mc phase near the triple point (grey dot at 
a=tan _1 (3), (3=tt) where the Mc, T, and O phases are 
in equilibrium. If a system such as BaTi03 or KNb03 
could somehow be perturbed so that the variation with 
temperature would carry the system on a trajectory pass- 
ing to the left of this triple point in Fig. 4, then a novel 
R-O-Ma-T-C (or even R-M B -0-M A -T-C) transition 
sequence might be observed. However, to our knowledge, 
no Mb or Mc phase has ever been observed in a cubic 
perovskite system. 



D. Higher-order expansions 

We have seen that the eighth-order expansion still does 
not allow for the appearance of a triclinic equilibrium 
phase for any parameter values. It is natural, then, to 
ask at what higher order in the expansion a triclinic phase 
can first occur. The answer is that the expansion must 
be carried to twelfth order before a triclinic phase can 
appear. 

For, suppose that a triclinic phase is the ground state. 
Then the minimum of G occurs at a point Uo located 
in the interior of the irreducible wedge, as illustrated in 
Fig. 3(d). Letting 

G(u) = G(u) - G(uo) , 

ff 4 (u) =3 4 (u)-3 4 (uo) , 



§6 (u) =5e(u)-5 6 (uo) , 
the expansion Eq. (|J) can be rewritten 

G(u) = 5 8 .g4 + 5io54.g6 + a 12 g\ + a' 12 gl + 



(12) 



(13) 



where the a n are trivial linear combinations of the a n . 
To drop the 0,4 and a^ terms, we have used the fact that 
G(u) must be stationary at Uo- (For this we also need 
that the gradients of g± and ge never vanish or become 
parallel at an interior point of the irreducible wedge; this 
is straightforward to confirm.) 

It is now evident that if expansion ( p~3[ ) is truncated at 
eighth order, then point Uo is not an isolated minimum. 
Instead, it belongs to a degenerate locus of minima cor- 
responding to 54 = 0, i.e., to a contour of the function 
<m(u). This is the situation illustrated in Fig. 3(e). It oc- 
curs when ae = 0, i.e., when cq + c§ = 0, corresponding 
to the grey vertical line at (3 — Zir/A in Fig. 4. 

If the expansion is carried to tenth order, then it is 
clear from Eq. ( ]13| ) that G still vanishes on this same 
contour. The two-dimensional Hessian matrix H^ = 
d 2 G / ' du^duy then takes the form 



H 



d 
d e 



(14) 



where indices fi=l and 2 correspond to the directions 
parallel and perpendicular to the 174 contour, respectively. 
d 7^ as long as aio 7^ 0, in which case det H < 0. 
Thus, at tenth order, the stationary point Uo cannot be 
a minimum; instead, it is generally an isolated saddle 
point. 

Finally, it is clear that the point Uo can be a local 
minimum if the expansion is carried to twelfth order. For 
suppose aio=5i2=0, and as and a' 12 are positive. Then 
G = asgl + ®i29e i s obviously positive definite, and Uo is 
a true isolated (global) minimum. 

Concluding this section, we find that it is necessary to 
go to surprisingly high order in the free-energy expansion 
in order to stabilize a ferroelectric state in which there 
are no symmetry constraints on the order parameter P. 
Specifically, we find that cubic invariants of twelfth or 
higher order have to be included to to stabilize such a 
triclinic phase. We conclude that the discovery (or syn- 
thesis) of a material having such behavior may be chal- 
lenging, but is by no means impossible. 



have undergone a symmetry-lowering structural distor- 
tion, but in such a way that long-range order has not set 
in. For example, the description of the phase transition 
sequence of BaTiOs and_KNb03 in terms of the well- 
known "eight-site model"E9 assumes the presence of ran- 
dom local rhombohedral displacements in the orthorhom- 
bic, tetragonal, and cubic phases. In the present case, 
Noheda et al. havft-concluded from their own structural 
analysis of PZTjLrti and that of Corker et al. of rhom- 
bohedral PZT,EU that there may be random local mon- 
oclinic displacements which order variously to yield the 
tetrahedral, rhombohedral, or monoclinic phases near the 
morphotropic phase boundary. One could then describe 
the system in terms of fluctuations between minima of a 
"24-site model." 

To make these ideas more precise, suppose that the 
local displacements of Noheda et al. arise from an optical 
branch of the phonon spectrum, and let u; be the vector 
"local made amplitude" for that branch within the Ith 
unit cell.EI Take as a model of the free energy 



F[ U ]=£/(uO + ^7zm(u,,U m ) 



(15) 



The in-cell energy / is presumed to be strongly nonlin- 
ear; the inter-cell coupling //,„ may either be bilinear in 
u/ and u m , or of higher order. / is also presumed to dom- 
inate the fi m so that, to a good approximation, its global 
minima establish the possible magnitude and orientations 
u" of u;. /(u/) can be subjected to precisely the same 
methods of analysis as applied to F(u) in Sees. ^-^V| 
above, with parallel results. The u" so obtained can 
then be substituted into F in Eq. (|15J), yielding. 



F[a] = F Q 



lm 



flmK, Ui) 



(16) 



This can be regarded as a 24-site version of the g-state 
Potts model, 

A statistical analysis can be carried out for vari- 
ous forms of fi m (e.g., bilinear) to capture the cu- 
bic (completely disordered), tetragonal (partially disor- 
dered), rhombohedral (partially disordered), and mono- 
clinic (fully ordered) phases of the case where u" takes 
on all of the 24 symmetrically equivalent monoclinic dis- 
placements. This type of analysis has already been car- 
ried out for orieatational order-disorder transitions and 
plastic crystals.Eij 



V. A MICROSCOPIC MODEL 

When structural transitions have some order-disorder 
character, a model free energy (e.g., (||) or (^j)) ex- 
pressed as a function of a macroscopic order parameter 
(e.g., u or u) provides little insight into the local struc- 
tural fluctuations that underlie the transitions. In such 
a case, a more appropriate microscopic picture of the 
high-symmetry phase may be one in which local regions 



VI. SUMMARY AND CONCLUSIONS 

The original Devonshire theorycl gave a natural expla- 
nation for the appearance of tetragonal, orthorhombic, 
and rhombohedral phases in materials such as BaTiOs 
based on a sixth-order free-energy expansion. Here, we 
have confirmed that these ferroelectric phases, in which 
the order parameter P is confined to a symmetry axis, 



are the only ones permitted by the sixth-order version 
of the theory. Moreover, we have clarified the nature 
of the phases that may be expected to appear at higher 
orders in the expansion. In particular, we have shown 
that the extension of the theory to eighth order allows 
one to describe, in addition, three kinds of monoclinic 
phase in which P is confined only to a symmetry plane. 
To obtain a triclinic phase in which P is unconstrained 
by symmetry, we have shown that a twelfth-order ver- 
sion of the theory is needed. A topological analysis of 
the critical points of the energy surface has been used to 
facilitate the discussion of the relevant phases and phase 
transitions. 

The present theory may provide some added insight 
into the phase behavior of conventional fcrroelectrics such 
as BaTi03, but the principal new results concern cases 
in which the eighth-order terms are important. In par- 
ticular, the theory provides a natural explanation for the 
monocluiiC|-|MA phase recently observed experimentally 
in PZTBEla It also predicts that if a triple point of equi- 
librium between T, R, and Ma phases occurs, then it will 
be one at which first-order R-T and R-Ma boundaries 
meet a second-order T-Ma boundary. That is, P will 
rotate into the (110) mirror plane continuously from the 
T side but discontinuously from the R side. 

Noheda et ala have shown such a triple point in their 
Fig. 6, in agreement with our analysis. However, they 
also show a region of coexistence of the T and Ma phases 
in the same phase diagram. In a homogeneous sample, 
such a coexistence region can be due to hysteresis aris- 
ing from nucleation barriers to a first-order phase transi- 
tion. The eighth-order theory predicts the T-Ma phase 
boundary to be of second order, in which case there 
can be no hysteresis or nucleation barriers. In princi- 
ple, higher-order contributions to the free energy could 
be large enough to change the order of the T-Ma tran- 
sition. However, a more likely explanation is that the 
samples studied by Noheda et al. may be inhomoge- 
neous. We note that these authors did not report hys- 
teresis but did report a two-phase coexistence region near 
the R-T boundary between 500 and 575 K in Fig. 5 of 
Ref. @ Attributing this to inhomogeneity, we estimate 
that the concentration variation may be of order 1% from 
the slope of the R-T/M-T phase boundary in Fig. 6 of 
Ref. fi Composition inhomogeneity of that magnitude 
would be sufficient to account for the T-MAphase coex- 
istence shown in that figure. 

Finally, the work may provide some guidance in the 
search for even more novel Mq, Mb, and triclinic ferro- 
electric phases. It also may be of utility in other kinds of 
cubic systems with other kinds of vector order parame- 
ters, e.g., ferromagnetic systems. 
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APPENDIX A: DETAILS OF SIXTH-ORDER 
THEORY 

The purpose of this Appendix is to give further details 
about the sequence of tra nsitions that occurs in the sixth- 
order model of Sec. IVB. Recall that the behavior in this 



model is governed by a single dimcnsionless parameter a 
defined via Eq. @. 

There are ten critical values of a that we can define as 

ai = -0.39767T = - tan _1 (3) 
a 2 =0.2256tt = tan" 1 (6/7) 
a 3 = 0.25tt 
a 4 = 0.3128tt = tan" 1 (3/2) 

«5 = 0.57T 

a 6 = 0.6024tt = 7r - tan" 1 (3) 
a 7 = 1.2256tt = 7T + tan" 1 (6/7) 
a§ = 1.25-7T 

a 9 = 1.3128tt = 7T + tan" 1 (3/2) 
otio = 1.57T (Al) 

The sequence of phases can be followed on Fig. 4 by trac- 
ing the vertical dotted lines, first from top to bottom at 
[3 = for < a < 7r, and then from bottom to top at 

(3 = 77 for 7T < a < 27T. 

In the T phase, the system exhibits the topology of 
Fig. 2(a) for u\ < a < a 2 ; Fig. 2(b) for a 2 < a < 
04; and Fig. 2(c) for a^ < a < a§. At a 2 , point O 
converts to a local maximum and simultaneously a saddle 
point appears along the Mb symmetry line. There is an 
irrelevant crossing of the R and O maxima at a^. At 
«4, R becomes a local minimum, and the saddle point 
switches from the Mb to the Ma symmetry line. The 
transition from T to R is first order at a§. 

The system falls into an R ground state corresponding 
to Fig. 2(c) for a§ < a < a%; Fig. 2(d) for a 6 < a < 07; 
and Fig. 2(e) for 07 < a < a 8 . The Ma saddle point 
disappears and T is converted to a maximum at a$, and 
a new Mb saddle point emerges with the conversion of O 
to a local minimum at 07. The transition from R to O 
is first order at as- 

Finally, an O phases occurs as illustrated in Fig. 2(e) 
for as < a < ag, and as in Fig. 2(f) for ag < a < ot\ +2n. 
The saddle point at Mb vanishes and R is converted to 
a local maximum at ag. There is an irrelevant crossing 
of the T and R maxima at aio- The transition from O 
to T at 1.6024tt is singular, in that the energy surface 
becomes exactly flat along the entire Mq symmetry line. 



1 A.F. Devonshire, Phil. Mag. 40, 1040 (1949). 
2 A.F. Devonshire, Phil. Mag. 42, 1065 (1951); Adv. Phys. 
3, 85 (1954). 

3 B. Jafle, W.R. Cook, and H. Jaffe, Piezoelectric Ceramics 
(Academic Press, London, 1971). 

4 M.J. Haun, E. Furman, S.J. Jang, and L.E. Cross, Ferro- 
electrics 99, 13 (1989). 

5 B. Noheda, D.E. Cox, G. Shirane, J.A. Gonzalo, L.E. 
Cross, and S.-E. Park, Appl. Phys. Lett. 74, 2059 (1999). 

6 R. Guo, L.E. Cross, S-E. Park, B. Noheda, D.E. Cox and 
G. Shirane, Phys. Rev. Lett. 84, 5423 (2000). 

7 B. Noheda, J.A. Gonzalo, L.E. Cross, R. Guo, S.-E. Park, 
D.E. Cox, and G. Shirane, Phys. Rev. B 61, 8687 (2000). 

8 B. Noheda, D.E. Cox, G. Shirane, R. Guo, B. Jones, 
and L.E. Cross, Phys. Rev. B, in press (preprint 



http://arXiv.org/abs/cond-mat/0006152D 



a W. Zhong, D. Vanderbilt, and K.M. Rabe, Phys. Rev. B 
52, 6301 (1995). 

10 L. Bellaiche, A. Garcia, and D. Vanderbilt, Phys. Rev. Lett. 
84, 5427 (2000). 

11 A.G. Souza Filho et al., Phys. Rev. B 61, 14283 (2000). 

12 Eq. (3) of Ref. hi] implies that the shear elastic constant 
vanishes at T = Tm even under conditions of fixed polar- 
ization P. Thus, the primary order parameter of the tran- 
sition is strain, not polarization, and the transition is to an 
"improper" or "extrinsic" ferroelectric state. 

13 For this to be so, the equilibrium magnitude |u| for a given 
u should not be too a strong function of u. 

14 F.C. Von der Lage and H.A. Bethe, Phys. Rev. 71, 612 
(1947). 

15 A. Garcia, private communication. 

16 M. Morse, Functional Topology and Abstract Variational 
Theory, Memoires Sciences Mathematiques (Gauthier- 
Villars, Paris, 1938), Fascicule 92; H. Seifert and W. 
Threlfall, Variationsrechnung Im Grossen, Hamburger 
Mathematische Einzelschriften (B.G. Teubner, Leipsig und 
Berlin, 1938), Vol. 24. 

17 J. C. Phillips, Phys. Rev. 104, 1263 (1956). 

18 L.D. Landau and E.M. Lifshitz, Statistical Physics, Part 1 
(3rd Edition) (Pergamon Press, New York, 1980). 

19 Representing u by polar coordinates (9, (f>) measured with 
respect to the (111) axis, a function such as sin 3 (#) cos(3<^>) 
is third order in displacements from 8 — and is invariant 
with respect to the rhombohedral symmetry group. 

20 R. Comes, M. Lambert, and A. Guinier, C.R. Acad. Sci. 
Paris 226, 959 (1968); Solid State Comm. 6, 715 (1968); 
J. P. Sokoloff, L. L. Chase, and D. Rytz, Phys. Rev. B 38, 
597 (1988). 

21 D.L. Corker, A.M. Glazer, R.W. Whatmore, A. Stallard, 
and F. Fauth, J. Phys. Cond. Matt. 10, 6251 (1998). 

22 W. Breymann and R.M. Pick, Europhys. Lett. 6, 227 
(1988). 



10 



